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Abstract 

The problem of glueball mass spectra using the relativistic Dirac equation is 
studied. Also the Breit-Fermi approach used to obtaining hyperfine splitting in 
glueballs. Our approach is based on the assumption, that the nature and the forces 
between two gluons are the short-range. We were to calculate the glueball masses 
with used screened potential. 
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1 Introduction 



The existence of gluon self-coupling in QCD suggests that, in addition to the conventional 
qq states, there may be non-qq mesons - bound states including gluons (glueballs and qqg 
hybrids). The existence of glueball states made from gluons is one of the important 
predictions of QCD. The discovery of these glueball states would be the strong support 
to the QCD theory. Therefore the search and identification of glueballs have been a very 
attractive research task. The abundance of qq mesons and possible mixing of glueballs 
and ordinary mesons make the current situation with the identification of glueball states 
rather complicated. So the spectrum of QCD is expected to contain glueballs hybrids 
and multiquark states. However, the theoretical quittance on the properties of unusual 
states is often contradictory, and models that agree in qq sector differ in their predictions 
about new states. Moreover, the abundance of qq meson states in the region 1 — 2GeV 
and glueball-quarkonium mixing makes the identification of the lightest non-qq mesons 
difficult. Whereas there is a general agreement that the lightest glueball should have 
quantum numbers ++ its mass remains controversial. The one of the possibilities to 
distinguish between possibilities is to study of production rates and decays branching 
ratios which are expected to be sensitive to the constituent structure. Glueballs are 
preferentially expected in production of gluon rich environment such as J/$ - , central 
production in hadron - hadron collisions by double Pomeron exchange and pp annihilation 
into mesons near threshold, or jj collision. 

The quark model predicts that there are two mesons with I = in the 3 Po qq nonet 
but apart from the states fj (1710) with J = or (and) 2, four states fo (400 — 1200), 
fo (980), / (1370) and fo (1500) are listed by Particle Data Group. There are too many 
controversies about these states. One of the possibilities is to put aside the f (400 — 1200) 
and fo (980) and focus on / (1370), / (1500), fj (1710) although the spin parity of the 
fj (1710) is J PC = ++ or (and) 2 ++ controversial. Recently several authors study the 
quarkonia-glueball contents of the /(1370), /(1500), fj (1710) by studying the mixing 
of these three states. The same situation is also in the case of identification of the 2 ++ 
glueball. Some progress has been made recently in the ++ scalar and 2 ++ tensor glueball 
sector, where both experimental and QCD lattice simulation results seem to converge. 

In this paper we will study the problem of glueball mass spectra using the relativistic 
Dirac equation, which was used to the calculations of mesons spectra, and the Breit-Fermi 
approach which was also used to obtaining the hyperfine splitting in different mesons. 
Our approach is based on the assumption that the nature and the forces between two 
gluons are the same that between two quarks, and because that the chosen approaches 
in both cases must be similar. In our previous papers we obtained a good description of 
meson spectra, fine and hyperfine splitting in relativistic Dirac approach with screened 
potential. Because that we were able to calculate the glueball masses and decays in the 
same approach, with the same parameters. 

2 Method Calculation of Glueballs states 

For calculation of glueball mass-spectrum we used Dirac equation with mixing Lorentz 
structure of potential for glueballs with spin averaged masses 
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- dJ j^ + $F ij (r) = [E-m-V}G ij (r) 
+ ±G i0 (r) = [E + m- V] F l3 (r) 



dr 

where 

7?7= ) =-(l + k). 

Let us write to equation all possible potential type 

d$> (r, t) 



= i^-V + (3 [m + V s (r) + 7 V P5 (r) + -fVf (r) + -f-fvf (r) + a m V^ 
The vector and scalar part we present in the following form: 



T 



V v (r) = Voce (r) + eV conf . (r) 
K s (r) = (l-e)V 0OR/ . (r) 



and substituted it to ([!]) 



dG(r) 
dr 



{E - Voce (r) - K on/ . (r) — m) F (r) = -f G (r) 
(£ - ^ogf (r) + (1 - 2e) K on/ . (r) + m) G (r) = ~F (r) + 



(E - Vqge (r) - V^/. (r) —m)F (r) 



dF(r) 
dr 



k 



r(E- Vqge (r) + (1 - 2e) K ori/ . (r) + m) 



(£ - V G gf (r) - Kon/. (r) —m)F (r) 



dF(r) 
dr 



r (£ - Voge (r) + (1 - 2e) Kon/. (r) + m) 
We obtain for the small component of the wave function G (r) : 

dG (r) 



dr 



(E - Voge (r) - Kon/. (r) —m)F (r) - 



r 



-£F(r) + 



dF(r) 
dr 



£ - Voge (r) + (1 - 2e) K on/ . (r) + m 



[(1 - 2e) K on/ . (r) - Vogf (r)] G (r) + 



dG (r) 

+ (E + (1- 2e) V conf . (r) - Vqge (r) + m) -^i = (9) 

= h-F (r) - -F' (r) + F" (r) . 

Into the formula (||) we substitute G (r), so we obtain the second order equation for the 
wave function F (r) : 

pa (r) _ l(l-2e)V conf .(r)-V OGE (r)]' 

y ' (E + (1 - 2e) V conf . (r) - Vfc GJS (r) +m) y ' 

, [ fc ^ , [(l-2e)V eon/ .(r)-y ocg (r)] / fc | 



r 



(£ + (1 - 2e) Kon/. (r) - ^ogs (r) + m) r 



+ (E + (1 - 2e) Kon/. (r) - ^oge (r) + m) (E - Vqge if) - V conf . (r) — m) } F (r) — 0. 

As for the interaction potential we choose the same screened potential as in the case 
of Breit-Fermi equation. With the same interaction potential we obtained good results in 
meson spectroscopy, and because that we try to apply this approach to glueballs. In the 
frame of Breit-Fermi approach the spin - spin interaction term is 

V ss = - — sls^AVy. (11) 

3m qi m q2 

We suggest for calculation to use configuration interaction approach, which was very 
successfully applied in atomic physics. 



3 Results 

In the Table 2 we show the results obtained in potential model for spin average glueball 
masses in comparison with the data of other authors. As it is seen from this table for 
the lowest glueball state our result is close for to the results of other authors. In our 
model the first and second orbital as well radial excited states are lower then in the 
lattice calculation |], 0. The same situation is also in the Table 3, where we show 
the results which incorporated the spin-spin effects. These probably mean, that the 
strong interaction between two gluons in glueball more strong than between two quarks 
in mesons. It is interesting to note, that to a same conclusion come J. Cui and H. Jin 0] 
during calculation of glueball masses, using Bethe-Salpeter equation. 



Table 1. Spin averaged glueball masses from Dirac equation with {V = 0.586GeV). 



Quantum 


Our results 


Nonperturb. method 


Lattice data 


Lattice data 


numbers 


M (GeV) 


[4] 


[1] 


[2] 


£ = 0,n r = 


4.68 


4.68 


4.66 


4.55 


£=l^n r = 


5.84 


6.0 


6.36 


6.1 


£ = 0,n r = l 


6.22 


6.0 


6.68 


6.45 


£ = 2,n r = 


6.69 


7.0 


9.0 


7.7 


£ = 1 , n r = 1 


6.96 


8.0 
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Table 2. Glueball masses from Breit-Fermi approach in comparison with other authors. 





our results, M (GeV) 


re/. [4] 


re/. [2] 


re/.[l] 


re/. [5] 


re/. [6] 


re/. [7] 


0++ 


1.70 


1.58 


1.73 


1.74 


1.645 


1.686 


1.659 


2++ 


2.19 


2.59 


2.40 


2.47 


2.337 


2.380 


2.304 


3++ 


3.24 


3.58 


3.69 


4.3 









The Decays of Glueballs to two light mesons: 
r (l 3 S 1 -> XX) = 55MeV, 
r (2 3 5i -> IX) = 24MeV, 
r (3 3 5i -> IX) = 15Me1/. 



4 Discussion and conclusion 



As it is note by Minkowski and Ochs ||, [|11|] in the gluon jet in case of triplet neutralization 
the leading particles are the hadrons formed by the primary gluon and soft qq pairs. 
If there are hybrid mesons they may be formed in the fast color singlet qqg system 
alternatively the fast hadrons are ordinary qq mesons formed at the end of the parton 
cascade after having absorbed all gluon energy. If the octet mechanism is at work the 
leading gluon may also form a glueball. On the other hand in the quark jet neither the 
hybrid nor the glueball will be leading if the leading quark is only neutralized in color by 
a soft antiquark. 

Table 3. Production of leading hadrons in the jet according H. 





neutralization 


qq 


hybrid 


glueball 


quark jet 


triplet 


yes 


no 


no 


gluon jet 


triplet 


yes 


yes 


no 




octet 


no 


no 


yes 



As it is note by Burakovsky || it is widely believed that pseudoscalar mesons are 
the Goldstone bosons of broken SU(3) xSU(3) chiral symmetry of QCD, and that they 
should be massless in the chirally symmetric phase. Because that it is not clear how the 
resonance spectrum would be suitable for the description of the pseudoscalar mesons. 

It is firmly established that the lightest glueball state is the scalar glueball. According 
to Burakovsky, who obtained the glueball masses from the Regge phenomenology, with 
linear Regge trajectory with negative intercept 



M 



T2 M{P) 



and 



M 2 (2 ++ ) = V2M (0 



So M (0++) = 1620 MeV, M (2++) = 2290 MeV . 

Tensor glueball is the lowest resonance lying on the Pomeron trajectory with unit 
intercept. As for this glueball state in PDG corresponds three candidates in this mass 
region fa (2220), J = 2 or 4, / 2 (2300) and / 2 (2340). 



4 



As to the question of hybrids. If hybrid spectroscopy works like conventional mesons 
spectroscopy or glueball spectroscopy than as shown by Toussaint |L(J, we expect that 



the lowest multiplet of hybrids contains 1~, 0~ + , 2 h particles both isovector and 
isoscalar which will be split by color hyperfme interaction. In this paper we restrict 
ourselves only with the calculations of glueball masses. 

Using a non-perturbative method based on asymptotic behavior of Wilson loops A. 
Kaidalov and Yu. Simonov calculated the masses of glueballs and corresponding Regge 
trajectories. These results are shown in the Tables 1, 2. 

The difficulties in identification of the glueball states, and the poor experimental data 
makes the investigation of the problem for glueballs actually. We try to obtain the mass 
spectrum of glueballs and the decay width in potential model approach. Using the same 
method and the same parameters as in the case of meson spectra calculation we obtained 
for glueball masses a good agreement with the possible experimental data. We obtained 
the results for spin average glueball states and the results which take into account the 
spin-spin interaction. Moreover, our results, as well as other authors results, shown, that 
the lowest scalar glueball mass is M (0 ++ ) = 1700 MeV, and the tensor glueball amass is 
M (2++) = 2290 MeV. Other results in Tables 1, 2 are predictions. 
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